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Ion Acceleration by Short Chirped Laser Pulses 
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Direct laser acceleration of ions by short frequency-chirped laser pulses is investigated theoreti¬ 
cally. We demonstrate that intense beams of ions with a kinetic energy broadening of about 1% can 
be generated. The chirping of the laser pulse allows the particles to gain kinetic energies of hundreds 
of Me Vs, which is required for hadron cancer therapy, from pulses of energies of the order of 100 J. 
It is shown that few-cycle chirped pulses can accelerate ions more efficiently than long ones, i.e. 
higher ion kinetic energies are reached with the same amount of total electromagnetic pulse energy. 
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I. INTRODUCTION 

The interaction of intense laser pulses with solids has recentl y a ttracted considerable interest. This is largely due 
to it^otential application for accelerating charged particles [lH22|. Tumor therapy with accelerated ion beams (see 
e.g. |23l-[^ j would in particular benefit from the replacement of conventional accelerators by all-optical devices, which 
may become compact and inexpensive in future, allowing this form of cancer therapy to be accessible for more patients. 

In this article, we demonstrate the feasibility of generating ion beams by shining an appropriately chirped short 
laser pulse on a target. At sufficiently high laser intensities the electrons are quickly ionized, and the ions get directly 
accelerated by the laser field. Modulating the frequency of the pulse leads to efficient particle energy gain from the 
field, as it was shown before [29143^. Here we consider the case of short chirped pulses, when the time duration of 
the pulse is comparable to a single cycle. We found that acceleration by such short pulses may be more efficient than 
by long chirped pulses, i.e. the same final ion kinetic energies can be reached with a lower pulse energy. 


Laser beam 

FIG. I: A schematic view of the studied laser acceleration set-up. A chirped few-cycle laser pulse is shot on a target. The 
ionized electrons are accelerated fast in the forward direction, which is followed by the more inert ions, also directly accelerated 
by the laser pulse. Note that the transverse target size is magnified for illustration purpose. In the simulations it is much 
smaller than the focal radius of the laser beam. 



The acceleration scheme is shown on Fig. [TJ We assume a solid-density or underdense target consisting of carbon 
ions, typically, bare nuclei, which can be generated via ionization by an intense pre-pulse. Collective plasma effects 
become important when the target thickness L in propagation direction of the laser is longer than the wave length 
of the induced plasma wave, i.e. Ap = 27rc/co’p with plasma frequency uj‘^ = Im^. In the last formula is the 

electron density, e the electron charge and the mass of an electron. When interacting with the accelerating pulse, 
first the lighter electrons are accelerated, they are pushed in the forward direction, as shown previously in Ref. [29|. 
The electrons are followed by the carbon ions, which are directly accelerated by the chirped pulse. 


II. LASER ACCELERATION SIMULATIONS 

In order to access the efficiency of laser acceleration by means of few-cycle chirped pulses, we perform simulations 
based on the classical relativistic equations of motion. In the following subsection, we provide a description of short 
laser pulses, based on earlier works [33|, [sj- Subsequently, we describe the particle dynamics in the presence of such 
intense chirped pulses. 
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A. Description of few-cycle chirped focused laser pulses 


We consider a circularly polarized laser field propagating along the ^-direction. In analogy with Ref. [35|, the 
vector potential for the focused ultra-short laser beam is represented as A = ^ (x'0(r, 77 ) + Here, 

the following notations have been introduced: the electric field amplitude is ^0 = 4>/P/c/7no, where P is the peak 
power of a laser pulse, c stands for the light velocity in vacuum, and wq is the laser focal radius. The phase of the 
field is given as 77 = ujot — k^z. The wave number is ko = ujq/c^ where ujq is the original (unchirped) frequency of 
laser at the focal point. When the chirping effect is taken into account, the laser frequency and radius have to be 
modified as follows: cj = cjq (1 + br]), and re = tuq/ (1 + 677 ), where b is the dimensionless chirp parameter. In these 
formulas, we introduced the quantities ^ = /(I + ir]/\ f = i/{i P vjZr\ z/ = 2 : + 77 / {2k), p = r/w, 
r = s = ujqt l2y/2 log 2 , and r is pulse duration. Furthermore, Zr = kw^ j 2 is the Rayleigh length, k = uj/c, 

and 00 is a constant phase. When introducing the frequency modulation, the phase parameter 77 and the pulse length 
parameter s do not need to be modified, x and y stand for unit vectors in the x and y directions orthogonal to the 
lasers propagation direction. 

Note that the temporal envelope of the laser beam is not factorized in this pulsed solution of the wave equation. 
The scalar potential 0 is thought to have a similar expression as that of the vector potential A, and it can be 
calculated from the Lorentz gauge condition 90/9t + V - A = 0. The electromagnetic fields are derived from E = 
—dAjdt — V0, B = \/ X A: E = E^^^ + E^y\ B = B^^^ + B^y\ herewith. 
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where Ei = V^o, with the normalization parameter So = [—si — l/{kzr)]/si, si = ^[1 + l/{2kzr) + l/s^], 

jz = z izr, and the following short-hand notations have been introduced: 


Cl 


C 2 


i + 


77 . 2^2 _ 2kjz — p + ips^ — 77 ^ 


+ ■ 


{2kjzEp)‘^ s^{p — is‘^) 

— 1 Akjz E 2p — 2ik‘^r‘^ p‘^ — s‘^ — ips‘^ 


{ 2 kjz + 77)2 


{ 2 k jz + 77)3 


{ips ■ 


^3\2 


2ip E s‘^ 
ips‘^ + 


( 3 ) 















































3 


These field expressions are in concordance with the sub-cycle pulse field of Ref. [36| and the long-pulse field of Ref. • 
Such fields were also employed recently to simulate quantum radiation reaction effects for an ensemble of electrons 
interacting with ultra-short pulses Q. For the much slower ions studied here, the influence of radiation reaction can 
be neglected [l5|. 


B. Particle dynamics 

The time-dependent dynamics of an ensemble of interacting ions is considered. An ion indexed by of mass m and 
charge q is accelerated to relativistic energy and momentum, respectively, of Ej = jjmc^ and pj = jjmcl3j^ where /3j 
is the velocity of the particle divided by c, the velocity of light in vacuum, and 7 = is the Lorentz factor, 

when interacting with the time-dependent Adds E and B of an laser pulse. Thus, the dynamics is described by the 
coupled Newton-Lorentz equations, given in SI (International System of Units): 

^ = q {E{rj) + • + c^j X {B{rj))) , (4) 

^ = qc(3r{E{rj) + E^;^ ) . 

The electric held of the inter-ionic interaction is approximated by with the Coulomb interaction 

scalar potential 


<t>jk = -^1——r- 

47reo |rj - Tfel 

Here, the relative displacement of two particles is Vjk = Vj — Vk and eo is the vacuum permittivity. The presence of 
plasma electrons can be neglected when simulating the ions’ acceleration dynamics, as they are blown off first by the 
pulse, as it was shown by particle-in-cell simulations [29|. 

Since the electromagnetic fields have a complex mathematical structure, and the inclusion of the ions’ interaction 
leads to coupled motion, one has to solve the above differential equations numerically. A numerical integration of 
Eqs. dl]) yields the particles’ velocities and thus also their energy gain at a Anal time equal to many laser cycles. 


III. RESULTS AND DISCUSSION 

We present first results of calculations for a single test ion. Bare carbon nuclei ions) are chosen here because of 
their higher contrast of dose deposition in the tissue as compared to protons. In order to And the optimal parameters 
of the pulse and the optimal initial velocity of the ion, we performed numerical simulations by calculating the energy 
gain for a particle as a function of the chirp parameter b and the initial energy Eq = jomc^ of the ion. The results are 
presented on Fig. [2l The optimal chirping is in the range around b = ±0.1, with low initial kinetic energies around 
7o ~ 1.05. We note that for the circularly polarized pulses employed here, the dependence of the energy on the 
chirping parameter is a slowly-varyin g fu nction. This is not the case for linearly polarized Aelds, where this function 
shows a strong oscillatory behaviour [^. Therefore, the appropriate chirping can be more practically implemented 
in experiments. The maximal energy gain which can be reached in this setting is around 100 MeV/u, reaching the 
range of interest in medical applications. 

As the previous Agure shows, the acceleration by a chirped short pulse is optimal when the initial kinetic energy 
of particles is low. Motivated by this, we separately consider the experimentally advantageous situation when the 
particles are initially at rest, i.e. 70 = 1. The energy gain is shown on Fig. [3] in dependence of the chirp parameter b. 
The maximal gain that can be achieved by such pulses when the particles are at rest initially is approx. 80 MeV/u. In 
the case of circularly polarized Aelds, similar gains can be reached when applying a negative frequency chirp (6 < 0 ), 
i.e. when the carrier frequency of the laser decreases in time. Additionally, circularly polarized laser beams have the 
advantage over linearly polarized beams that they are not sensitive to small variations of the chirp parameter b - cf. 
Fig. [3] with Fig. 2 of Ref. [29|. Fig. 0] shows, for an optimal chirp parameter b = 0.089, the energy gain and the 
corresponding pulse as a function of the longitudinal displacement z, conArming that the Anal kinetic energy around 
100 MeV is reached on a sub-wavelength scale. The Agure also illustrates - in accordance with earlier Andings [29| 
- that it is the asymmetric part of the pulse, induced by the strong frequency modulation, which causes sudden 
acceleration of the ion. This is in contrast to an unchirped laser pulse, where a charged particle would gain kinetic 
energy during the Arst half of a laser cycle and subsequently loses this energy during the second half of the respective 
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FIG. 2: Energy gain in units of MeV/u in dependence on the chirp parameter b and the initial dimensionless energy parameter 
7o (i.e. Lorentz factor, Eq = yomc^). The wavelength of the laser is chosen to be 1 /rm and its power to 10 PW. The pulse 
duration is equal to 3 laser cycles, and the focal waist radius is equal to A. A circularly polarized laser pulse is employed. 



FIG. 3: Energy gain per nucleon for G®^ nuclei as a function of the chirp parameter b. The nuclei are assumed to be initially 
at rest, i.e. 70 = 1. The other parameters are as described for Fig. [2 


cycle due to the symmetry of the pulse. However, we want to emphasize that frequency modulation over a large 
frequency bandwidth is presently experimentally accessible only at lower field intensities [^, . 

In order to find the optimal pulse duration for acceleration, we performed single-particle simulations calculating the 
energy gain of a particle initially at rest in dependence of this parameter, at a fixed total pulse energy. The results 
are displayed on Fig. [5] (a), showing that maximal efficiency can be reached for pulses with a time duration of 3 to 4 
cycles. Increasing the pulse duration can only decrease the final energy gain, since the peak electric field of the pulse 
decreases if the energy of the pulse is distributed over a longer time. Furthermore, in pulses with durations below 3-4 
periods, the ion does not spend sufficient time interacting with the field to reach its maximal velocity. For each pulse 
duration, the chirp parameter b was optimized independently. The corresponding values of the chirp parameters are 
displayed on Fig[5](b). 

The energy distribution of a beam of ions generated by laser acceleration was determined by many-particle calcula¬ 
tions based on the coupled Newton-Lorentz equations [Eq. dH]. We assume that the ions are randomly distributed in 
a small nano-scale cylindrical volume before interacting with the short focused laser pulse. The ions were assumed to 
be initially at rest, and have a solid-state density (2.4x10^^ ions/cm^). Simulations were performed with and without 
the inclusion of the ion-ion interaction to assess the effect of Coulomb repulsion on the energy gain and its broad¬ 
ening. In both cases, the set of random initial coordinates utilized was kept constant. Results are shown on Fig. [6l 
For this case collective plasma effects are negligible since the induced plasma wavelength, Xp = 0.0088/im, is much 
larger than the target thickness, L = 0.0004/im. Including the particle-particle repulsion yields an order-of-magnitude 
broader distribution than the simulation with artificially switching off the interaction, however, the average kinetic 
energy is not influenced by this: in the first, realistic simulation, the average gain is 82.557 MeV/u, with a standard 
deviation of 1.042 MeV/u (1.262%); in the latter case, the average gain is the same, while the standard deviation is 
0.006 MeV/u (0.007%). Thus, we may conclude that for dense targets, it is indeed necessary to include inter-particle 
interactions in realistic simulations, even at the very high laser intensities considered. The beam energy spread is in 
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FIG. 4: (a) Energy gain per nucleon for nuclei and (b) the electric field of the pulse as a function of the phase 77 , for a 

fixed optimal chirp parameter of 5 = 0.089. The other parameters are as described for Fig. O 



FIG. 5: (a) Energy gain per nucleon for G®^ nuclei and (b) and optimal chirp parameter 6 as a function of the pulse duration 

T, given in units of the period tq. In these calculations, the pulse energy was kept fixed at a value of 100 J. 


the medically applicable range, i.e. approximately 1 %, enabled by the small size of the target. Trajectory of the 
overdense interacting ion beam is illustrated in Fig. [71 and the emission angle is approximately 0.05 (radian), i.e. 
2.865 (degree). Please note that the angle employs the unit of radian in the rest content. For the same conditions 
except taking a larger chirp parameter, b = 0.1, the emission angle increases from 0.05 to 0.051. And, for a longer 
laser pulse with the pulse duration r = 4ro, the emission angle decreases to 0.048. 

Besides a nano-scale target, simulations were made for a more extended underdense carbon plasma target, which 
may be realized as an expanding plasma created by a pre-pulse. The kinetic energy distribution for a target of ions 
with an initial density of 2.4x10^^ ions/cm^ is presented as a histogram on Fig. |8l The induced plasma wavelength, 
Xp = 0.88/rm, and the target thickness, L = 0.04/im. Therefore, collective plasma effects are ignored. In this case, 
particle-particle interactions are weaker, and thus the results are less sensitive to the presence of interaction. We 
obtain an average ion gain of 82.553 MeV/u, with a standard deviation of 0.138 MeV/u (0.167%) for the more 
realistic simulation with the ions’ mutual repulsion taken into account, and to the same average gain with a spread 
of 0.069 MeV/u (0.083%) with the interaction neglected. The energies reached here coincide with the results for a 
solid-density ion plasma target, however, the energy broadening is somewhat better. Trajectory of the underdense 
interacting ion beam is illustrated in Fig. |9l Emission angle is approximately 0.005 and 10 time smaller than 
that of the overdense case. 
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FIG. 6 : (a) Energy distribution histogram for an ensemble of 3000 interacting ions, confined initially in a small cylinder 

(length 0.0004 /xm, diameter: 0.002 /xm). (b) Energy distribution for the same set of ions, with the inter-particle Coulomb 
repulsion switched off. Note the different energy scales of the two sub-figures. 
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EIG. 7: Trajectory of the overdense interacting ion beam. Red curves are the trajectories of ions, and blue points 
are the corresponding initial and final positions in the simulations with — 4s < 77 < 4s. The parameters employed here are the 
same as those of Eig. [ 6 ] (a). 


IV. CONCLUSIONS 

Relativistic many-particle simulations were performed in order to access the applicability of few-cycle chirped laser 
pulses to ion acceleration for the purpose of hadron cancer therapy. An accurate description of the temporal and 
spatial structure of the laser fields was employed. We have found that such pulses with durations of 3-4 cycles, when 
focused on small underdense or solid-density plasma targets, can produce ion beams with properties in the range of 
medical requirements. They are also more efficient in accelerating ions than longer chirped pulses investigated in our 
earlier work [29|: a further increase of pulse duration beyond the optimal value decreases the particles’ energy gain, 
as the pulse energy distributed over longer times leads to weaker accelerating fields. 
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FIG. 8: (a) Energy distribution histogram for the case of an underdense target formed by 3000 interacting ions, conhned 

initially in a cylinder (length 0.04 /xm, diameter: 0.02/im) (b) Energy distribution for the same underdense target, with the 
inter-particle interaction neglected. 



EIG. 9: Trajectory of the underdense interacting G®^ ion beam. Red curves are the trajectories of G®^ ions, and blue points 
are the corresponding initial and hnal positions in the simulations with —4s < rj < 12s. The parameters employed here are the 
same as those of Eig. [8] (a). 
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